Kondo Screening Cloud and Charge Quantization in Mesoscopic Devices 
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We propose that the finite size of the Kondo screening cloud, £x, can be probed by measuring 
the charge quantization in a one-dimensional system coupled to a small quantum dot. When the 
chemical potential in the system is varied at zero temperature, one should observe charge steps that 
are controlled by the Kondo effect when the system size L is comparable to £k- We show that, if 
the standard Kondo model is used, the ratio between the widths of the Coulomb blockade valleys 
with odd or even number of electrons is a universal scaling function of £,k /L. 

PACS numbers: 72.15.Qm, 73.21. La, 73.23.Hk 



The Kondo effect can be described as the strong renor- 
malization of the exchange coupling between an electron 
gas and a localized spin at low energies [lj|. Below some 
characteristic energy scale, known as the Kondo temper- 
ature Tk, a non-trivial many-body state arises in which 
the localized spin forms a singlet with one conduction 
electron. In the recent realizations of the Kondo effect Q , 
a semiconductor quantum dot in the Coulomb blockade 
regime plays the role of a spin S = 1/2 impurity when 
the number of electrons in the dot is odd. The Kondo 
temperature in these systems is a function of the coupling 
between the dot and the leads and can be conveniently 
controlled by gate voltages. As a clear signature of the 
Kondo effect, one observes that the conductance through 
a quantum dot with S = 1/2 is a universal scaling func- 
tion G (T/Tk) and reaches the unitary limit G = 2e 2 /h 
when T <C Tk for symmetric coupling to the leads [3|. 
On the other hand, we expect that in a finite system the 
infrared singularities that give rise to the Kondo effect are 
cut off not by temperature, but by the level spacing A if 
A ^> T. For a one-dimensional (ID) system with length 
L and characteristic velocity v, the relevant dimension- 
less parameter is A/Tk ~ €k/L, where £k = v/Tk is 
identified with the size of the Kondo screening cloud, the 
mesoscopic sized wave function of the electron that sur- 
rounds and screens the localized spin. This large length 
scale (£k ~ 0.1 — l^m) is comparable with the size of cur- 
rently studied mesoscopic devices, which has motivated 
several proposals that the Kondo cloud should manifest 
itself through finite size effects in such systems [3, Q . 

One familiar property of small metallic islands in the 
Coulomb blockade regime is the quantization of charge. 
Since each electron added to the system costs a finite 
energy, the number of electrons changes by steps as 
one varies the chemical potential and sharp conductance 
peaks are observed at the charge degeneracy points [H- 
In AlGaAs-GaAs heterostructures, the level spacing of a 
ID system with size L ~ 1/j.m is of order A ~ lQQ/ieV, 
large enough to be resolved experimentally. Clearly, the 
energy levels should be affected by the Kondo interaction 



FIG. 1: Possible experimental setup. Vdw controls the tun- 
neling t' between the small dot (on the left) and the wire and 
V g varies the chemical potential in the wire. 



with an adjacent dot. Therefore, the addition spectrum 
should exhibit signatures of the finite size of the Kondo 
cloud. 

We consider the geometry shown in Fig. [TJ One of the 
ends of a wire of length L with a large number of elec- 
trons (such that A <C £f, where cf is the Fermi energy) is 
weakly coupled to a small quantum dot via a gate voltage 
Vdw ■ Both the dot and the wire (which can be thought of 
as a very long and thin dot) are in the Coulomb blockade 
regime. The dot has a very large charging energy and 
a ground state with S = 1/2. We assume that the wire 
is very weakly connected to a reservoir and that trans- 
fer of a single electron between wire and reservoir could 
be measured by some Coulomb blockade technique. This 
coupling is assumed weak enough so as not to affect the 
energy levels of the wire-dot system. Experiments with 
small dots and large dots have been performed, for ex- 
ample, in Q and considered theoretically in [f|. 

Here we investigate this effect using the basic Kondo 
model, ignoring charge fluctuations on the dot, assum- 
ing the wire contains a single channel and ignoring all 
electron-electron interactions in the wire. Within the 
constant interaction approximation these interactions 
are represented by a simple charging energy whose effects 
on the charge steps are easily subtracted off. We post- 
pone a more complete treatment of these complications, 
to a later work. (The effect of short range interactions on 
the conductance through a Kondo impurity was treated 
in the Luttinger liquid framework in [5(.) 
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Weak coupling. - The tunneling between the wire and 
the dot is well described by the Anderson model [B| . For 
zero tunneling amplitude the system reduces to an open 
chain and a free spin. In the continuum limit, we in- 
troduce the fermionic field ^(x) for the electrons in the 
wire, as well as its right- and left-moving components 

* ( x ) = e tkFX ^ R (x) + e~ ikFX ip L (x) . 

The open boundary conditions (0) = (L) = im- 
ply that tpL,R are not independent. Instead, ipL can be 
regarded as the extension of tpR on the negative- a; axis: 
ipR (—x) = —ipL (x). We can then work with right movers 
only and write down an effective Kondo model (we drop 
the index R hereafter) 

pL -. 

H=-iv F dxtp f d x ^ + 2nv F X^ (0) (0) • S, (1) 
J-l 2 

where A is the dimensionless Kondo coupling. We assume 
that A is independent of V g , which may require that t' 
and the energy of the localized state in the dot be tuned 
accordingly. For A <C 1, the size of the Kondo screening 
cloud (defined in the thermodynamic limit) is exponen- 
tially large: £k ~ (vp/D) e 1 ^, where D -c e F is a high- 
energy cutoff. The boundary condition in the weak cou- 
pling fixed point (L/£k — > 0) reads vp(—L) = e lS if;(L), 
where 5 = 2kpL mod tt. Using the mode expansion 
ip (x) = — («'/v2i) J2k e%kx<z ki we obtain the momentum 
eigenvalues k n = (nir — S) / L, n e Z. 

We denote by N(fi) the total number of electrons in the 
system, including the one in the quantum dot. We are in- 
terested in the elementary steps of N around some initial 
value N = N (jug). At T = 0, we calculate N (fj,) as the 
integer value of N that minimizes the thermodynamic po- 
tential Q (TV) = E (TV) - fiN, where E (TV) is the ground 
state energy. Defining as the value of ijl where TV 

changes from Nq + £ to No + £ + 1, it follows from the 
charge degeneracy condition (N + £ + 1) = ft (N a + £) 
that fj, e+ x = E {N + £ + 1) - E (N + £). For A = 0, 
we set fio halfway between the highest occupied and the 
lowest unoccupied energy level of the open chain. In 
this case, iVo = odd: the ground state is doubly de- 
generate (total spin S tot = 1/2) and consists of one 
electron in the dot and pairs of conduction electrons in 
the Fermi sea. The energy levels measured from ^o are 
e„ = (n— 1/2) A, where A = nvp/L. It is easy to see 
that, as we raise /x > N jumps whenever (i crosses 
an energy level. Moreover, there are only double steps 
due to the spin degeneracy of the single-particle states 



(/ i 2m+i — A^m+f)- A s a result, N (fi) ~ No takes on 
only even values and TV is always odd (see dashed line in 

Fig. rj. 

We obtain the charge steps in the weak coupling limit 
by calculating the correction to the ground state energy 
using perturbation theory in A Q. In the limit A — > 0, 
the ground state for TV odd is still doubly degenerate and 
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FIG. 2: (Color online) Charge quantization steps for the wire 
coupled to a small quantum dot. The arrows indicate the 
direction in which the single steps move as A e / / (L) grows. 



we write \GS (N = odd)) = n„< 4141 1°) ® It) , where 
It) = i i s the spin state of the dot. However, for 
N even, there is one single electron occupying the state 
at the Fermi surface (with momentum hp). To lowest 
order in degenerate perturbation theory, the ground state 
for A is \GS{N = even)) = Il„<o4I c feI l°) ® l s )> 
where |s) ee [\k F |) ® \fy) - \k F {) ® |fr)]/V2 is the singlet 
state between the spin of the dot and the electron at 
kp. This state has Stot = and ( S ■ S e i) = —3/4, 



where S e 
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c k is the total spin of the conduction 



electrons. The singlet formation lowers E (N) for N = 
even relatively to N = odd. As a result, the Kondo 
interaction splits the double steps and gives rise to small 
plateaus with N even (Fig. [2J. To O (A 2 ), we find 
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Vp 



(2) 



where [£] = £ + 1 for £ even and [£] = £ for £ odd. Note the 
logarithmic divergence at O (A 2 ) as L ~ * oo, character- 
istic of the Kondo effect. We recognize the expansion of 
the effective coupling A e // (L) ~ [In (^/L)] -1 in powers 
of the bare A, as expected from scaling arguments. We 
have verified that the scaling holds up to 0(A 3 ), i.e., 
Mf+i has no dependence on the cutoff D but the implicit 
one in the expansion of A e / f (L). This is remarkable given 
that E (TV) itself is cutoff dependent. Based on this, we 
conjecture that the relative width of the plateaus with 
TV odd/even is a universal scaling function of £,k/L. We 
define the ratio (see Fig. [5]) 

_ 6fj° _ E(Nq + 3) - 2E(N + 2) + E(N Q + 1) 

~ 6fi e E(N Q + 2) - 2E(N + 1) + E(N ) ' ( ' 

From Eq. @, we have that in the weak coupling limit 



3 V L / 3 \ L 



-const.. (4) 
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Strong coupling. - When X e ff — * oo, the spin of the 
dot forms a singlet with one conduction electron and de- 
couples from the other electrons in the chain. In this 
limit the Kondo cloud is small (£k/L — ► 0). The strong 
coupling boundary conditions reflect the tt/2 phase shift 
for the particle-hole symmetric case: ip (-L) = —ip (L) 
Q. This implies that the k eigenvalues are shifted with 
respect to weak coupling: k n = (mr + Tr/2 — 5)/L. At this 
fixed point we recover the spin degeneracy of the single- 
particle states and double steps in N (fi) appear at the 
shifted energy levels (e„ = raA as measured from the 
original /Xq ) . The ground state consists now of a singlet 
plus pairs of electrons in the wire, thus N (fi) is always 
even. Fig. [2] illustrates how the charge staircase evolves 
monotonically from weak to strong coupling. 

We explore the limit fcj? 1 C ^ < I by working out 
a local Fermi liquid theory 0]. The idea is that virtual 
transitions of the singlet at x = induce a local interac- 
tion in the spin sector of the conduction electrons. The 
leading irrelevant operator that perturbs the strong cou- 
pling fixed point and respects SU(2) symmetry is 



H 



FL 



3 T K 



(O)-ifr(O) 



(5) 



where the prefactor is fixed such that the impurity sus- 
ceptibility is Ximp = 1/(47a'). This interaction lowers 
the ground state energy when there is an odd number of 
remaining conduction electrons (S e i = 1/2), thus split- 
ting the charge steps in the strong coupling limit. To 
lowest order in I/Tr, we find 



A 2 8 L '' 



from which we obtain 
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Bethe Ansatz results. - In order to check the scal- 
ing property of R = 5 ^i odd / 5 fi cvcn for any L/£k we 
use the Bethe Ansatz (BA) solution of the one-channel 
Kondo problem We start with a half-filled band 

of N — 1 conduction electrons (N odd) coupled to a 
localized impurity spin, corresponding to a system size 
L = (No — l)/2 in units where the Bethe ansatz cut- 
off parameter D [l(|, related to the band- width, is set 
equal to one. Since the BA solution can be obtained 
for any filling factor, we can add particles one by one 
in the system (TV = No, Nq + 1, N + 2, . . .) and com- 
pute the corresponding energies. This has been done by 
solving numerically the coupled BA equations (Io| . us- 
ing a standard Newton-Rapson method. In Fig. [3] we 
present results obtained for the ratio R = 5 /i odd / 5 /j, even 
with iVo = 51, 101, 201, 501, 1001, 2001 and 13 different 
values of the Kondo coupling in the range 0.06 < J < 0.8. 
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FIG. 3: (Color online) Universal ratio R = S^ odd /6^ cvcn as a 
scaling function of £k/L. Bethe ansatz results obtained for 
various systems sizes (N = 51, 101, 201, 501, 1001, 2001) 
and 13 different values of the Kondo exchange J, indicated by 
different symbols. For each value of J, the system lengths L 
have been rescaled L — > L/£k(J) in order to obtain the best 
collapse of the data using the strong coupling curve Eq. J7} 
(dashed red line) as a support for the rest of the collapse. 
The weak coupling regime for L <C £x, enlarged in the in- 
set (a), is described by the weak coupling expansion Eq. Q 
with constant ~ 0.33 (continuous blue curve). Inset (b): The 
Kondo length scale, extracted from the universal data col- 
lapse of the main panel (black squares), is described by the 
exponential fit Eq. ((SJ) with £o — 0.31 (dashed green line). 



The universal scaling curve has been obtained by rcscal- 
ing the x-axis L — > L/£k(J) in order to get the best 
collapse of the data. The entire crossover curve is ob- 
tained, ranging from the strong coupling L £k to the 
weak coupling regime L <ti £k- Both strong coupling 
[Eq. [7J] and weak coupling [Eq. results are perfectly 
reproduced. As displayed in the inset (a) of Fig. O BA 
results agree with Eq. (U) using only one fitting param- 
eter: constants 0.33. The Kondo length scale, shown in 
the inset (b) of Fig. [3l displays the expected exponential 
behavior II Oil 



Tr/C 



(8) 



with c = 2,7/(1 — 3/4J 2 ). Fitting £r- to this expression 
gives £o — 0.31 (see inset (b) of Fig. [3]) which is in very 
good agreement with the expected value of £o = I /it, 
resulting [l[ from the impurity susceptibility normalized 
to 1/(4T K ). 

Flux dependence in quantum rings. ~ The effect we 
have described is not unique to the wire geometry. One 
interesting alternative is to suppose that the quantum 
dot is embedded in a ring with circumference L (inset of 
Fig. rjj • This geometry offers the possibility of looking at 
the dependence of the charge steps on the magnetic flux 
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FIG. 4: (Color online) Flux dependence of the charge steps 
in the ring with an embedded quantum dot (inset), assuming 
zero potential scattering. Dashed lines: weak coupling regime 
(£k 2> L); solid lines: strong coupling regime (£jf -C £)• 



4> threading the ring, which is related to the persistent 
current for the embedded quantum dot [B|, [ijj . For sim- 
plicity, we neglect potential scattering and assume that 
the coupling between the dot and the leads is symmetric. 
In this case it is convenient to label the eigenstate of the 
open chain by their symmetry under the parity transfor- 
mation x — > L — x. We define the even and odd fields 
^ ( x ) = ( x ) ± e~ lkFL ^ L (L - x)]/V2, in terms of 
which the Kondo interaction takes the form Q 



a 



Ot X , . (X ,X , 

cos-tpl (0) -isva-ipl (0) 



Ot Qt 

cos-V e (0)+isin-V o (0) -5, (9) 



where a = 27r$/$ and = hc/e is the flux quantum. 
For Nq = 4p+ 1 (p integer), the lowest unoccupied state 
of the open chain (A = 0) is even under parity. Defining 
5fx n = M2n+2 — M2n+i as the splitting of the charge steps 
in the weak coupling limit, we find 

^ = \ [1 + (-1)" cos a] X eff (L) + O (\j ff ) . (10) 

Notice that for a = 2m-K only the even channel couples 
to the impurity and the steps in the odd channel (corre- 
sponding to odd values of n) do not split. The opposite 
happens when a = (2m + l)ir (see Fig. |4j. 

In the strong coupling limit, the 7r/2 phase shift modi- 
fies the boundary conditions in such a way that the trans- 
mission across the dot becomes perfect 0, [j| . The system 
is then equivalent to an ideal ID ring, whose eigenstates 
can be labeled as right or left movers. The energy lev- 
els are flux dependent and exhibit zigzag lines with level 
crossings at a = mir 12|. The major effect of a Fermi 
liquid interaction analogous to Eq. © is to lift the de- 
generacy between right and left movers at a = mir and 




to split two out of the four values of /ig+i [13|. We calcu- 
late the splitting 5/j,' n = ^2n+^ ~M2n-§ using degenerate 
perturbation theory to 0(1/Tr-) and And (£k <C L) 



(11) 



for a sa m7r,m = n,n ± 2, . . . (Fig. 2J. This striking 
difference in the flux dependence of the charge steps is 
consistent with the results in 

Finally, we would like to point out that the constant 
interaction model assumed here neglects the exchange in- 
teraction between electrons in the wire. However, such 
interaction could be cast into one empirical parameter 
that gives a constant contribution to the even-odd split- 
ting, even in the weak and strong coupling limits 0j. 

In conclusion, we have shown that the width of the 
charge steps in a mesoscopic device depends on the finite 
size of the Kondo cloud. The signature of this effect 
would be the broadening of the Coulomb blockade valleys 
with total N even as we increase Tjc (decrease The 
crossover between weak ^> L) and strong 
L) coupling could be used to extract the characteristic 
length scale of the Kondo problem. 
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